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Abstract. We present the light curve of an unusual vari- 
able object, DUO 2, detected during the search for mi- 
crolensing events by the DUO project. The star remained 



qv^ stable for more than 150 days before it brightened by more 
. than two magnitudes in 6 days in the B and R bands. The 
^2 light curves are achromatic during the variability. We con- 
H-j sider possible explanations of the photometric behavior, 
(*-■>) with particular emphasis on the binary lens interpretation 
£S) of the event. The masses of the lenses are quite small, with 
the companion possibly in the range of a brown dwarf or 
even a few times of Jupiter. We report evidence of blending 
of the source by a companion through the first detection 
^ — I of shift in the light centroid among all the microlensing ex- 
periments. This shift sets a lower limit of 0.3" on the sep- 



aration between the stars. The best lens model obtained 



requires moderate blending, which was what motivated us 
to check the centroid shift that was subsequently found. 
The best lens model predicts a separation of 1" between 
I the two blended stars. This prediction was recently tested 
O using two CCD images taken under good seeing condi- 
> tions. Both images show two components. Their separa- 
tion and position angle are in good agreement with our 
model. 



spectacular case is the first binary lensing candidate dis- 
covered by Udalski et al. (1994b), which was later con- 
firmed by Alcock et al. (1995). Lensing by binary systems 
is naturally expected in about 10 percent of the lensing 
events (Mao & Paczyhski 1991). 

The DUO (Disk Unseen Objects) collaboration is a 
project with the main goal of searching for galactic dark 
matter (Alard et al. 1995). One of the microlensing candi- 
dates discovered presents a very unusual variability. In this 
paper, we present a detailed analysis of this event, with 
particular emphasis on the binary lensing interpretation, 
much in the same spirit of Udalski et al. (1994b). 

2. OBSERVATIONS 

2.1. data collection 

Data presented here were collected during the 1994 ob- 
serving season with the ESO lm Schmidt telescope. 

The Schmidt telescope provides photographic plates 
of 28cm x 28cm. The scale is 67"/mm, giving a wide field 
of 5.2x5.2 degrees. Two B plates (IllaJ) and one R plate 
(IIIaF) were taken for each photometric night. The mean 
time interval between two B plates was 1.5 hours, with 
the R plate taken just in between. 



Key words: dark matter - gravitational lensing - stars: 
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1. INTRODUCTION 

The search for dark matter using gravitational microlens- 
ing (Paczyhski 1986) has been impressively successful, 
with several groups discovering many microlensing can- 
didates (Bennett et al. 1994, Alcock et al. 1993; Aubourg 
et al. 1993; Udalski et al. 1994a; Alard et al. 1995). One 
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2.2. reduction 

All the images were reduced with the standard package 
of the DUO project (Alard 1995). The analysis of light 
curves revealed an unusual microlensing candidate. This 
event displays a multi-peaked structure while maintain- 
ing the achromaticity. Table 1 contains the coordinates, 
magnitudes, and color for this event. Because the star is 
close to one corner of the plate, we emphasize that the 
magnitude zero point could have been affected. We also 
caution that the colors are seriously affected by the high 
and variable extinction in the region surrounding the star. 
The 60 x 60 arcsecs region centered on the star is shown in 
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Figure 1 (Plate ESO11307B). This can serve as a finding 
chart. 



Table 1. Positions and magnitudes for DUO 2 
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Fig. 1. 60 arcsec by 60 arcsec B-band image centered on the 
variable star. Two black bars identify the lensed star. North is 
up and east is to the left. 



3. Nature of Variability 

We now explore the possible explanations for the multi- 
peaked variability exhibited by this object, including vari- 
able stars and lensing by binaries. 

3.1. An Eruptive or Flare Star? 

The rapid rises and falls in the light curve suggest some 
type of eruptive or flare stars (cf. Udalski et al. 1994b for 
a similar discussion). The flare star possibility seems very 
unlikely as all the known flare last no more than a few 
hours and there is no other flaring activity in our data 
which may be expected for a flare star (cf. Udalski et al. 
1994b). A dwarf nova or a low-mass X-ray binary (LMXB) 
can also explain the observed amplitude of observation. 



However, we note that the candidate is close to the bulge 
turnoff point where most stars are located, i.e., there is 
nothing unusual in its B-R color. While the cataclysmic 
variables usually have some UV excess, the U-B color for 
this star does not differ from the colors of the majority 
of the stars. Furthermore, the flat baselines in B and R 
are unusual for a cataclysmic variable. We also searched 
for the outbursts of LMXBs in the IAU circulars for the 
observational period, none was found in the direction of 
the source. 

Future optical spectral observations will help to re- 
veal the nature of the object, as dwarf novae usually have 
strong emission lines in the minimum state. Furthermore, 
if the variability repeats itself in the future, the eruptive 
nature of the object will be established. 

3.2. A Double Lens? 

The first binary event discovered in the direction of the 
bulge (Udalski et al. 1994b; Alcock et al. 1995) shows that 
binary lensing events can indeed be identified in the mi- 
crolensing experiments, as predicted by Mao & Paczyriski 
(1991). In comparison with the first event, our light curve 
has one more distinctive peak. Our "daily" time coverage 
translates to a worse coverage relative to the duration of 
variability. Nevertheless, the achromaticity in the B and R 
bands and the unusual variability both suggest that this 
may be due to a binary lens, a possibility we shall explore 
below. 

For our binary lens model, there are 11 parameters: 
the total mass of the binary, M; the mass ratio of the in- 
dividual masses, q = m,ijm,\\ their separation, a, as pro- 
jected onto the lens plane; the closest approach to the 
center of mass, b; the time of closest approach, tf,] the an- 
gle between the axis of the binary and the trajectory, 9; 
the magnitudes of the star at the minimum light in the 
B and R bands, nriB and nriR, and finally the finite size 
of the lensed source, r s . Two more parameters, fn and 
Jb , specify the light contributions at the baselines of the 
lensed source in R and B. We assume the source surface 
brightness can be described by a limb darkening profile, 
I\(r) oc 1 — u\ + u\y^l — r 2 /r^, where r s is the radius 
of the star, and u\ is a wavelength dependent parameter 
(Allen 1973). We adopt ub = 0.6 and ur = 0.5, the values 
for a solar type star. The source profile has little influence 
on the results as it is not well constrained due to sparse 
sampling. 

The fitting of a binary light curve has been studied 
by Mao & Di Stefano (1995). The fitting first finds good 
initial parameters and then uses these good guesses as in- 
put to search for acceptable fit(s) by minimizing x 2 ■ F° r 
our candidate, we searched for the good initial guesses 
by Monte Carlo simulations. For this, we fix the baseline 
magnitudes at tur = 18.55, tub = 20.21, as these pa- 
rameters are well determined. Furthermore, we ignore the 
source size. The mass ratio, separation, impact parame- 
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ter, and angle are generated uniformly in the range of: 
q = 0.l- 1, a/r E = 0.5 - 2.5, b/r E = - 1.25, = - 2tt. 
For each set of q, a, b, and 0, the light curve was computed. 
Those events with at least three peaks are further consid- 
ered. For these events, we adjust the two time parameters, 
t E and t\,, such that two of the multiple peaks in the the- 
oretical light curves coincide with the first two observed 
ones. The initial parameters with x 2 < 348 for 116 data 
points are selected for further minimization. For this final 
step, we vary all the parameters to search for acceptable 
fits. 

In Fig. 2 we plot our best fit model. The model has a 
X 2 of 89 for 116 data points (105 degrees of freedom). The 
best fit parameters are: 

q = 0.33; a/r E = 1.21; b/r E = 0.40; = 94°.62; 

t E = r E /V = 8.5 day; t b = 85.4 day; f B = 0.73; f R = 0.70; 

r s /r E <, 1.0 x 10" 2 ; m B = 20.21; m R = 18.55; (1) 

The model reproduces the observed light curve very well. 
The lensing geometry is shown in Fig. 3. The trajectory 
has two caustic crossings and passes close to two cusps, 
producing four peaks in total. The first peak is missed due 
to sparse sampling. 



Fig. 2. The B and R band light curves for DUO 2. The solid 
line shows the light curve of the best binary microlens model. 
Two insets show the light curves of the lensed parts in two 
bands. The parameters for this fit are described in eq. (1). 



3.3. Blending 

As we showed in the last section, the best lensing model re- 
quires a moderate (30%) light contribution of an unlensed 
component. Originally, we only explored models without 
blending. For such models, we found consistent 3er — 5er 
deviations, especially for the point around day 85, where 



Fig. 3. The geometry of the best fit binary microlens model. 
The two components are labeled as two black dots. The caustics 
(thick solid line) and critical curves (dashed line) are shown. 
The caustics are the source positions where the magnification 
formally diverges for a point source. The caustics on the source 
plane are mapped into critical curves in the lens plane. The 
trajectory of the source is indicated as a straight line with an 
arrow indicating the direction of the motion. Both axes are in 
units of the Einstein radius (te ~ 1 AU) in the source plane. 



the magnification is 3.1 ± 0.15. If we assume the two ad- 
jacent peaks are produced by two caustic crossings, then 
such low magnification requires rather special parameters 
to achieve - as the absolute minimum magnification is 3 
when the source is inside the caustics (Witt & Mao 1995). 

The blending requirement in the model prompted us 
to perform an additional check on the images. If the ex- 
tra light is provided by a chance superimposition in the 
crowded field, then one may expect a positional shift in 
the centroid of light. The centroid position can be deter- 
mined typically to 0.2 pixel (0.66"/pixel) for B and 0.35 
pixel for R at the minimum state. For the bright state, the 
determination improves to about 0.1 pixel and 0.2 pixel, 
respectively. The shift of the centroid, Ar c , from that of 
the minimum state, should follow an explicit relation: 

Ar c = (l-/)r {I- A- 1 ), (2) 

where, r is the vector from the position of the unlensed 
component to the position of the lensed component. In 
Fig. 4, we plot the positional shifts in units of pixel (0.66") 
in the B band versus (1 — A -1 ). By linear regression, we 
obtain (cf. eq. [2]). 

(1 - f)Ax c , B = -0.41 ± 0.09, (1 - f)Ay C}B = 0.21 ± 0.09; 

(l-f)Ax c , R = -0.42±0.23, (l-f)Ay c , R = 0.10±0.23,(3) 

where x increases from east to west, and y from north to 
south. The correlation is statistically significant in the B 
band but somewhat marginal in R, although the two are 
consistent. As 1 — / < 1, eq. (3) provides a lower limit, 
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0.3", on the separation between the two blended compo- 
nents. The orientation of the two components is 60° west 
of north. For the best lens model, / « 0.7, a separation of 
1" is then inferred. 

After the completion of our analysis, two CCD I band 
images of the object were taken under seeing of 0.8" (Van 
Der Hooft 1995; Szymahski & Udalski 1995). The compo- 
nents are resolved due to the good quality of the images. 
The location of the unlensed star is in very good agree- 
ment with our model. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



1 — 1 /A 1 - ' /A 

Fig. 4. The B band shifts in the x and y direction of the 
centroid of light are shown versus the 1 — A^ 1 , where A is the 
linear increase in flux (magnification). The shifts are expressed 
in units of pixel (0.66"). The solid lines are obtained by linear 
regression (eq. [3]). 



4. DISCUSSION 

In the previous sections we have presented the peculiar 
light curves of DUO 2. The binary lens model can well 
reproduce the observed light curve. The model requires 
moderate blending, which was later confirmed by the shift 
in the light centroid, and the CCD images. The confir- 
mation of our prediction on the position of the blending 
component with the the CCD images does increase our 
own confidence in the lensing interpretation. The achro- 
maticity of the light curve further supports this interpre- 
tation. Notice that although the two components are likely 
to have somewhat different colors, as the unlensed com- 
ponent is only contributing a moderate fraction of light, 
the achromaticity is still roughly maintained. Due to the 
short time duration of the variability and our short obser- 
vation season, we cannot rule out the possibility that this 
star is just an unusual variable star, i.e., the lens model 
is simply a false positive (cf. Mao & Di Stefano 1995). If 
this event turns out to be a variable, the importance of 
frequent sampling is then observationally highlighted. Al- 
though we have searched a reasonable range of parameter 
space, we can not guarantee that the fit we presented is a 
unique one. 

Let us now assume that the binary lens interpretation 
is correct and proceed to estimate the mass of the lens and 



the size of the source. The total mass of the lens is given 
by (e.g., Paczynski 1986) 

m = »- 013m °GooiM ! (t^) (^) ! ' (4) 

where x = Dl/Ds is the ratio of the distances to the lens 
and the lensed star, Vj is the transverse velocity projected 
onto the source plane and we have adopted Ds =8kpc. 
The transverse velocity and the distance to the lens are 
both uncertain, but it seems likely that the lenses are low 
mass objects and the less massive lens can easily have 
the mass of a brown dwarf or even just a few times that of 
Jupiter. Such low mass lenses are unlikely to contribute to 
the observed light. The size of the lensed source is not well 
constrained, as the peak behaviors are not well resolved, 
only a rough upper limit can be obtained: 

Such a limit can be easily satisfied for a star like our source 
near the turn-off point. 

Regardless of whether the variability is induced by a 
binary lens, we demonstrated, for the first time in all the 
microlensing experiments, that blending can be detected 
through carefully examining the shift in the centroid po- 
sition. 
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